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Abstract. We consider a sample of 22 nearby clusters of 
galaxies observed with the Medium Energy Concentrator Spec- 
trometer (MECS) on board BeppoSAX . They cover the range 
in gas temperature between 3 and 10 keV, with bolometric X- 
ray luminosity between 2 x 10 44 erg s _1 and 6 x 10 45 erg s _1 . 
Using the de-projected gas temperature and density profiles re- 
solved in a number of bins between 5 and 7 and obtained from 
this dataset only, we recover the total gravitating mass profiles 
for 20 objects just applying the (i) spherical symmetry and (ii) 
. . . hydrostatic equilibrium assumptions. We investigate the corre- 
lations between total mass, gas temperature and luminosity at 
' several overdensities values and find that the slopes of these re- 
: lations are independent of the considered overdensity and con- 
sistent with what is predicted from the cluster scaling laws. The 
best-fit results on the normalization of the M — T relation are 
slightly lower, but still consistent considering the large errors 
^-j- that we measure, with hydrodynamical simulations. A segre- 
t-H gation between relaxed and non-relaxed systems is present in 
each plane of these relations pointing out a significant compo- 
nent in their intrinsic scatter. This segregation becomes more 
, evident at higher overdensities and when physical quantities, 
' like M gas and L, that are direct functions of the amount of gas 
y— ( , observed, are considered. 

\o ■ 

Key words: galaxies: cluster: general - galaxies: fundamental 
' parameters - intergalactic medium - X-ray: galaxies - cosmol- 
ogy: observations - dark matter. 
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1. Introduction 

The amplitude and the shape of the power spectrum of the pri- 
mordial density fluctuations on scales of about 20 Mpc 
can be effectively constrained with the mass function of galaxy 
clusters. Since the early '90s, X-ray observations have been 
used to build large datasets of measured luminosities and, with 
more effort because a larger number of source counts is re- 
quired, temperatures of the X-ray emitting plasma trapped in 
the cluster gravitational potential. These observed quantities 
are expressions of the physical processes that are taking place 
in the galaxy clusters and manifest the energy and the mass 
of these systems. Then, comparing the observed distribution in 
luminosity (or temperature) with theoretical models of the ex- 
pected cluster number density that are functions of total mass 
and redshift and depend upon the cosmological model adopted 
(e.g. Press & Schechter 1974), it has been possible to put con- 
straints in the "normalization-shape" plane of the primordial 
density fluctuations spectrum (see, e.g., the pioneering work of 
Henry & Arnaud 1991 and the most recent results in Ikebe et 
al. 2002 and references therein). 

However, the conclusions reached making this comparison 
rely on an efficient way to relate the observed quantities (like 
gas luminosity and temperature) to the gravitating mass of the 
systems. Gas-dynamics simulations (e.g. Evrard, Metzler & 
Navarro 1996, Schindler 1996) have confirmed the expected 
correlation between mass and temperature and have shown that 
mass estimates are reliable when obtained through X-ray anal- 
ysis under the assumption of spherical symmetry and hydro- 
static equilibrium. More recently, mass profiles obtained relax- 
ing the condition of plasma isothermality have shown a sig- 
nificant mismatch in normalization and slope of the mass- 
temperature relation between observational data and simula- 
tions (e.g. Horner, Mushotzky & Scharf 1999, Nevalainen, 
Markevitch & Forman 2000). On the other hand, it has been 
clear since the first compilation of catalogues of luminosity and 
temperature (Mushotzky 1984, Edge & Stewart 1991) that the 
observed correlation between these two quantities deviates sig- 
nificantly from the expected scaling law, suggesting contribu- 

Send offprint requests to: S. Ettori 
Correspondence to: settori@eso.org 



tions to the total energy of the plasma from physical phenom- 
ena other than the gravitational collapse. 

The observed Luminosity-Temperature (L — T) and Mass- 
Temperature (M — T) relations for galaxy clusters are, there- 
fore, the foundation to construct the cluster mass function and 
to use these virialized objects as cosmological probes. In this 
paper, we investigate these relations and, more in general, 
any correlation between observed and inferred quantities us- 
ing BeppoSAX observations of 22 nearby clusters of galaxies 
with resolved gas temperature and density profiles. The main 
differences between this study and previous work on the same 
subject are: 

1 . the use of BeppoSAX data that allows us to extend the anal- 
ysis of spatially-resolved spectra up to 20' in radius, i.e. ~ 
2.5 times the most favourable configuration with Chandra 
(Weisskopf et al. 2000) and to put under control some sys- 
tematic effects (e.g., sharper and more energy-independent 
Point-Spread-Function than the ASCA one -Tanaka et al. 
1994-, more stable and lower background than the one ob- 
served in XMM-Newton -Jansen et al. 2001), 

2. the direct deprojection of the spectral results to recon- 
struct the gas temperature and density profiles in a model- 
independent way. 

The sample presented in this work is, to date, the largest for 
which the physical quantities (i.e. gas density, temperature, lu- 
minosity, total mass, etc.) have all been derived simultaneously 
from spatially-resolved spectroscopy of the same dataset. The 
difference between this approach and others which make use 
of data coming from different satellites and/or make strong as- 
sumptions on the temperature profiles, such as isothermality, is 
twofold: on one side the use of data from different missions and 
the simplistic assumptions on the temperature profiles allow to 
build up samples bigger than ours, on the other they increase 
the likelihood of systematic effects which may in turn affect 
the relations between the observed quantities. 

The paper is organized as follows: in Section 2, we describe 
the BeppoSAX MECS observations of the galaxy clusters in our 
sample and the results of the spectral analysis considered in this 
work; the deprojection technique applied to the projected spec- 
tral results is discussed in Sect. 3; the gravitating mass profiles 



S. Ettori et al.: Cluster mass profiles using BeppoSAX 



3 



are obtained and compared with the optical measurements in 
Section 4; in Section 5, we study the correlation between the 
total mass, gas temperature, gas mass and luminosity; we sum- 
marize our results and present our conclusions in Sect. 6. 

All the errors quoted are at la level (68.3 per cent level 
of confidence for one interesting parameter) unless otherwise 
stated. The cosmological parameters Ha = 50/i^~ km s _1 
Mpc -1 and £l m = 1 — Cl\ = 1 are assumed hereafter. 



2. The sample 

We selected from the BeppoSAX SDC archive all the on-axis 
pointings of galaxy clusters with redshift smaller than <~ 0.1 
and exposure times larger than 30 ksec. The observation log 
for the cluster sample, with a detailed discussion of the metal 
abundance and temperature profiles derived for subsets of this 
sample, is given in De Grandi & Molendi (2001 and 2002, here- 
after DGM02). 

In this paper we discuss data from the imaging Medium- 
Energy Concentrator Spectrometer (MECS; 2-10 keV; Boella 
et al. 1997). The MECS consists of two identical grazing in- 
cidence telescopes with imaging gas scintillation proportional 
counters in their focal planes. The field of view of the MECS 
is circular with radius of <~ 25'. This detector has a spectral 
resolution of <~ 8% at 6 keV and a Point Spread Function 
(PSF) of ~ 1' (HPR), which varies only weakly with the en- 
ergy (D'Acri, De Grandi, & Molendi 1998). The MECS has 
an entrance Beryllium window sustained by a thicker support- 
ing structure, or strongback, in form of a circular ring and four 
ribs, which has transmission properties different by the rest of 
the window. 

The data analysis is fully described in De Grandi & 
Molendi (2001, 2002), hence in this paper we will only sum- 
marize the whole procedure. Standard reduction procedures 
and screening criteria have been applied using the SAXDAS 
package under the FTOOLS environment to produce equal- 
ized and linearized MECS event files. Each cluster has been 
divided into concentric annuli centered on the X-ray emission 
peak computed by fitting a Gaussian to the photon distributions 
in both the x- and y-direction on ROSAT PSPC images; out to 
8' we accumulate spectra from four annular regions each 2' 
wide; beyond this radius we accumulate spectra from annuli 4' 
wide. The energy dependent PSF of the MECS and the energy- 
dependent telescope vignetting for on-axis observations have 
been taken into account in our extended sources analysis by 
generating appropriate instrument response files (with the ef- 
farea program available within the SAXDAS package) to be 
used when fitting the accumulated spectra. We have computed 
the corrected effective area for the 8' — 12' annulus, which is 
covered by the circular region of the strongback, by consider- 
ing the typical thickness of the strongback and its transmission 
as a function of the energy and position. All other regions of 
the detector covered by the strongback have been appropriately 
masked and the data rejected. The background subtraction has 
been performed using spectra extracted from blank sky events 
files in the same regions of the detectors as the source. 



We fitted each spectrum with a single-temperature plasma 
in collisional equilibrium at the redshift of the cluster (Mekal 
model -Kaastra 1992, Liedhal et al. 1995- in XSPEC v. 10.0 
-Arnaud 1996), absorbed by the nominal Galactic column den- 
sity (wabs model; Dickey & Lockman 1990). 

3. Deprojection of the spectral results 

The physical quantities constrained from fits of spectra with 
counts collected from cluster regions projected on the sky need 
to be converted to their values in the spherical shells that con- 
stitute the assumed spherical geometry of the X-ray emitting 
plasma. Fitting a thermal model to a projected spectrum pro- 
vides, for each annulus, (i) an estimate for the Emission Inte- 
gral, EI = J n c n p dV = 0.82 / n\dV , through the normal- 



ization K of the model, K 



icr la 



EI (see Mekal 



model in XSPEC; we assume n p = 0.82n o in the ionized 
intra-cluster plasma); (ii) a direct measurement of the emission- 
weighted gas temperature, T r ; ng (note that the observed temper- 
ature is properly a photon-weighted temperature, but the dif- 
ference from our assumed definition is completely negligible), 
metal abundance, Z r [ ng , and luminosity, L r i ng . The purpose of 
the deprojection is, for example, to recover the value of the gas 
temperature in shells, T s h n = T if that is defined as 



T • = T- — 



outer shell 



TiWi 



(i) 



outer shell W K 



where Wy = Li x Vol(i,j)/Vol(i) = CiVol(i,j) provides 
the luminosity for a given shell i with volume Vol(i) weighted 
by the part of this volume projected on the ring j, Vol(i,j). 
Using this notation, it is simple to note that L r i ng = Lj = 

Si, outer shell e i^ol(i,j) — outer shell Wi 3 • 

From Kriss, Cioffi & Canizares (1983; see also McLaugh- 
lin 1999 and, particularly relevant to X-ray analysis, Buote 
2000), the volume shell observed through each ring adopted 
in the spectral analysis can be evaluated and a matrix, Vol, can 
be built with components equal to the parts the volume of the 
shells (rows i) seen at each ring (or annuli; column j). 

The deprojected physical quantities can be then obtained 
through the following matrix products (shown by the symbol 
#): 



;/, = (Vol T )- 1 #(S//0.82) 

e = (Vol T )- 1 #L riM 
eTkeii = (Vol T )- 



1/2 



4^ ( ^ r i n g T r i n g 



(2) 



eZshell = (Vol ) 1 #(L r ing^ring), 



where (Vol T ) _1 indicates that the matrix is firstly transposed 
and then inverted. The emission due to the shells projected 
along the line of sight but with the corresponding annuli out- 
side the field-of-view is taken into account with an edge correc- 
tion factor estimated assuming a power law distribution of the 
emission proportional to r~ 4 (cfr. equation A8 in the appendix 
of McLaughlin 1999). Finally, we have to assign a single ra- 
dius, r avo , to each shell. Formally, for each shell delimited from 
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the radii r, and r^+i, this radius should be the one that solves 
the equation n 2 gas {r avc ) = 3/(rf +1 - rf) f^ +1 n 2 gas (r)r 2 dr. 
Considering that (i) n gas (r) cx r~ a with a that has generally 
a value enclosed between 3 and 6 and can also vary between 
these values as function of radius in the same cluster, and (ii) 
complicated iterative procedure and interpolation are required 
(see, e.g., discussion about equation A9 in the appendix of 
McLaughlin 1999), we have checked that the assumption of 
^avo = (fi+t + is consistent with analytic solution of the 
equation above for an acceptable a within 2 per cent. 

We have applied this technique to the single-phase results 
of the spectral analysis presented in DGM02. In rings where 
part of the flux was masked for the presence of point-sources, 
we correct the normalization K by the relative amount of area 
not considered implicitly assuming spherical symmetry. 

For each cluster, we finally have the following outputs: gas 
bolometric luminosity in each shell, L = £ s hcii; gas tempera- 
ture in each shell, T = T s h n; electron density in each shell, 
n e , and, integrating it over the volume, the gas mass, M gas . An 
error is assigned to each quantity given the distribution of the 
values after 100 Monte-Carlo simulations obtained from scat- 
tering the original projected input with respect to their Gaus- 
sian error. 

Following DGM02, we divide our sample into two groups 
of objects with (CF; 12 clusters) and without (NCF; 10 clus- 
ters) a cooling flow in the central region (see reviews in Sarazin 
1988, Fabian 1994), according to the mass deposition rate 
quoted in Peres et al. (1998; NCF systems have a deposition 
rate consistent with zero; cf. Table [j]). It is worth noticing, how- 
ever, that this classification is meant to distinguish between 
relaxed (CF) and not-completely- relaxed (NCF) systems. In 
fact, also considering recent results from Chandra (David et 
al. 2001, Ettori et al. 2002, Sanders & Fabian 2002, Johnstone 
et al. 2002) and XMM-Newton (Tamura et al. 2001, Molendi & 
Pizzolato 2001, Matsushita et al. 2002) analyses of nearby clus- 
ters that show a lack of spectroscopic evidence of multi-phase 
gas, we parameterize all the X-ray emission with a single tem- 
perature model. 

4. X-ray estimate of the gravitational mass profiles 

To estimate the total gravitating mass, M to t, we make direct 
use of the deprojected gas temperature and electron density 
values estimated from the spectral best-fit with a single phase 
model. For each cluster, we select the mass model that repro- 
duces better the deprojected gas temperature profile inverting 
the equation of the hydrostatic equilibrium between the dark 
matter potential and the intracluster plasma: 
n e M to t,modci(< r) _d(n c x kT) 



Gfj,m p - 



(3) 



r 2 dr 
where /i=0.6 is the mean molecular weight in a.m.u., G is the 

gravitational constant, m p is the proton mass, and using the 
deprojected electron density, n c . As mass models, we consider 
two functional forms obtained from the integration of the fol- 
lowing dark matter density profile: (i) the King approxima- 
tion to the isothermal sphere (King 1962, Binney & Tremaine 
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Fig. 1. Histograms of the number of objects considered at each 
overdensity A in our analysis. The solid line shows the number 
of clusters where the outer radius, i? ut, to which the quantities 
in exam are observed is larger than (or equal to) r& at the 95 
per cent level of confidence [i.e. i? ou t > (?"a — 1.96c), where 
a is the error quoted in Table ^ as function of the overdensity 
A. The dotted line indicates the total number of clusters in our 
sample. The dashed line shows the number of objects that do 
not satisfy the selection criterion in mass (i.e. M to t > at any 
given r a at the 95 per cent level of confidence; see Sect. 5). 
The number of objects considered at each A is given, therefore, 
from (dotted line - dashed line). 



1987), with a flat core in the inner part and ar 3 dependence at 
r — > oo; (ii) the function discussed in Navarro, Frenk & White 



(1997, hereafter NFW), with a r" 1 and a 
in the inner and outer parts, respectively: 

M to t,modcl(< r) = 47T p s f(x), 



r 



-2.4 



200 



= Pc 



3 ln(l + c)-c/(l + c) 



dependence 



(4) 



m = 



ln(x + Vl + x 2 ) 



-n^f (King) 
(NFW) 



Hi + x)-^ 

where x — r/r s , p c is the critical density and the relation 
?"A=200 — c x r s holds for the NFW profile. 

Both of these mass models have two free parameters, the 
core (King) or scale (NFW) r s and the normalization, that we 
quote through the concentration parameter c (note that we do 
this also for a King profile for convenience). For each mass 
model, we obtain the best-fit values of the two parameters min- 
imizing the x 2 °f the comparison between the deprojected tem- 
perature profile and the one obtained from eqn. |] in two suc- 
cessive steps. First, a minimum in a x 2 distribution is searched 
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Table 1. Sample of 22 galaxy clusters considered in this study (see DGM02 for details on the reduction and spectral analysis of 
the BeppoSAX observations). The 'CF' column indicates if a cluster is or not a cooling-flow system according to the amount of 
deposition rate quoted in Peres et al. (1998). i? ut is the value of the radius at the outer end of the last radial bin considered. The 
best-fit parameters r s (in kpc) and c for a King and a NFW dark matter density profile are quoted with the respective errors (in 
parentheses). Note that A1367 and A3376 do not converge in the r s — c plane. 
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varying these parameters within the following ranges: 10 kpc 
< r s < max(2000kpc, -R ou t), 1 < c < 15. This search provides 
the best-fit values r' s , c'. A second fit is, then, performed in the 
restricted ranges: [min(rg-300 kpc, 10 kpc), max(rg+300 kpc, 
Rout)], [min(c'-3.0, 0.5), c'+3.0]. The results of this refined fit 
on the scale radius r s and the concentration parameter c are 
presented in Table |l|. Hereafter, M to t(< f) is defined for each 
cluster according to the minimum \ 2 provided from the two 
mass models considered, Af to t,King(< t) and M to t,NFw(< r )- 
The error related to the mass estimate is obtained from half the 
difference between the maximum and the minimum value cal- 
culated at each radius for the set of parameters acceptable at 1 
a. 

From our final sample, we exclude A1367 and A3376 be- 
cause we do not obtain any \ 2 solution for them. Out of the 
remaining 20 objects (12 CF and 8 NCF systems), ten (6 of 
which are CF clusters) are fitted better with a King profile. 

We investigate the relations among different physical quan- 
tities considering their values at a given overdensity, A. 
This is defined with respect to the critical density, p c z = 
(iHl)/ (8irG), and within a cluster described as a sphere with 
radius r&: 

A = 3M t ot(< r A ) 



with the Hubble constant at redshift z equal to 

H z = H Q Vfi m (l + zf + 1 - fi m (6) 

(for Q m + fl\ = 1; i.e., Hq X (1 + z) 3 / 2 for an Einstein-de 
Sitter universe). 

The following analysis has been performed at different 
overdensities. To handle the observed profiles at any radius, 
we interpolate linearly all the quantities on scales of 1 kpc. 
In Fig. [j], we show the number of clusters for which the re- 
gion enclosing a given overdensity A is directly accessible to 
our X-ray observations. From this figure, we conclude that, at 
A =2500, 18 galaxy clusters have a detectable X-ray emis- 
sion and two (A426 and A3526) need an extrapolation of the 
physical quantities (RA/Rout — 1-19 and 1.34 for A426 and 
A3526, respectively). As reference value for our results (cf. 
Table Q), we consider also A = 1000, where 1 1 objects are 
observable, eight (A85, A119, A426, A496, A1656, A2029, 
A3526, A3571) need an extrapolation in radius by about 40 
per cent (5, 29, 87, 15, 39, 11, 95 and 20 per cent, respec- 
tively) and one (A3627) did not satisfy our selection criterion 
in mass (a^j/M =0.62 > 0.51, see Sect. 5). Furthermore, to 
compare our results with previous work, we estimate the quan- 
tities examined, i.e. gas density, temperature and luminosity, 
at lower overdensity. For those clusters without observed val- 
ues at these ta, we extrapolate the interesting quantities us- 
ing a least squares error-weighted fit with a first-order polyno- 
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Table 2. Results from the deprojection analysis. All the quantities are estimated within J?a (apart from T, which is estimated at 
Ra), where the given overdensity A is obtained assuming either a King or NFW functional form for the total mass profile (la 
error in parentheses). 
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4.05 (0.43) 


6.01 (0.24) 
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2.17 (0.30) 


5.25 (0.54) 
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A426 
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1078 (80) 
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9.30 (0.82) 
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1321 (36) 


3.69 (0.44) 


5.73 (0.14) 


5.15 (0.25) 


A2029 


1561 (32) 
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7.67 (0.28) 
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8.03 (0.38) 
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4.43 (0.12) 


4.41 (0.21) 
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5.71 (0.30) 
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8.96 (0.73) 
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2.73 (0.13) 
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7.51 (0.18) 
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6.67 (0.85) 
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19.37 (0.49) 
6.61 (0.50) 
37.61 (0.39) 
6.54 (0.19) 
23.96 (0.82) 
22.20(1.19) 
23.91 (0.44) 
46.63 (1.12) 
58.18 (1.26) 
7.46 (0.18) 
19.12(0.72) 
38.44(1.87) 
20.62(1.14) 
2.67 (0.12) 
4.24 (0.39) 
19.10(0.53) 
10.29 (1.42) 
7.61 (0.18) 
65.11 (0.92) 
36.81 (1.45) 



11.33 (0.25) 
6.97 (0.67) 
9.09 (0.07) 
4.24 (0.08) 
15.23 (0.17) 
10.61 (0.16) 
10.11 (0.19) 
17.91 (0.27) 
19.89 (0.38) 
4.09 (0.13) 
13.50 (0.22) 
17.89 (0.32) 
13.01 (0.51) 
1.78 (0.05) 
3.86 (0.36) 
13.18(0.23) 
4.77 (0.29) 
3.47 (0.10) 
15.81 (0.21) 
15.17(0.49) 



9.14 (0.61) 
10.44 (1.88) 
4.74 (0.27) 
2.90 (0.26) 
10.32 (0.57) 
3.90 (0.87) 
8.05 (0.65) 
13.80 (0.86) 
10.90 (0.80) 
2.87 (0.19) 
9.65 (0.95) 
10.16(1.83) 
9.98 (0.93) 

I. 00 (0.13) 
2.76 (0.59) 

II. 73 (1.21) 
2.76(1.70) 
2.27 (0.16) 
11.78 (0.88) 
9.10(1.32) 



mial performed with the svdf it function (Press et al. 1992, 
sect. 15.4) on the logarithmic values of the variables observed 
in the outer region where r > 0.7 x i? ou t- The mean relative 
error measured in the observed region is propagated to the ex- 
trapolated values. 

For a given overdensity A, we quote in Table ^ the values 
of Af to t(< ?"a), ?~a, several estimates of the gas temperature 
(see Sect. 5), L bo i(< r A ) and A/ gas (< r A ). 



4.1. X-ray mass: comparison with (3"/-model 

In this section, we compare the estimates of the dark matter 
profile we have obtained in the previous subsection with results 
derived from modelling the gas density profile and applying (i) 
the hydrostatic equilibrium and (ii) a polytropic shape of the 
temperature profile. The latter procedure is generally applied 
in the X-ray analysis of galaxy clusters and makes use of the 
j3— model, p gas oc (1 + x 2 )~ 1 ' 5 ^ (Cavaliere & Fusco-Femiano 
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1976), to reproduce the observed surface brightness profile (an 
analytic expression can be obtained if T gas is assumed constant; 
for a generalization to the polytropic case see Ettori 2000) and 
build a temperature profile as function of the polytropic index 
7 like T gas oc (1 + ^-i.s^-i). 

In the /37-model the total mass profile is readily derived 
from eqn. |3| 

Af tot ,/9 7 (< r) 0.6 x 3 

-- 1.11 — p 7 T(r) r c 



Mass (<r 1000 ) 



10 14 h^M ( 



{l + x 2 ) 



1.11 — P 7 lb r c 



(7) 

This formula has been applied to estimate the total mass 
profile in recent work that considered a measured tempera- 
ture profile from ASCA data (e.g. Markevitch et al. 1999 on 
A496 and A2199; Nevalainen, Markevitch & Forman 2000, 
Finoguenov, Reiprich & Bohringer 2001). 

In fig. ||, we compare our mass profiles (from the King func- 
tional form, cf. Table [l]) for two CF clusters (A496 and A2199) 
with those derived by Markevitch et al. (1999) using the j3j 
model. The larger deviations (|om| > 3er) are localized in the 
region ~ 100 — 500 kpc (and below 100 kpc and above 1 Mpc 
in A496) and introduce a systematic error that could contribute 
to the observed scatter in the distribution of the measurements. 
It is worth noticing that the polytropic temperature profile does 
not reproduce in a satisfactory way the temperature profile of 
either A496 or A2199, as can be seen in fig. §. DGM02 have 
shown that the temperature profiles of our BeppoSAX sample 
of clusters are in general not in good agreement with poly- 
tropic temperature profile. This discrepancy should be taken 
into consideration when applying the f3j model to derive mass 
measurements of galaxy clusters. 

4.2. X-ray mass: comparison with optical estimates 

Girardi et al. (1998) quote the optically-determined mass es- 
timates for 15 out of 20 of our clusters (not available for 
A3526, A3627, 2A0335, PKS0745, TRIANG). In Fig. § we 
show a comparison between X-ray measurements at A =1000 
and the optical masses estimated at the same r& by using the 
Jeans equation (Mj SO = S^^a^r / G, where /3 ga i is the ex- 
ponent in the King-like galaxy density profile and a p is the 
projected velocity dispersion. This mass estimate is consis- 
tent with the corrected virial mass as discussed in Girardi et 
al. (1998, Section 5). Moreover, by making use of the relation 
a 2 = (GM200)/ (2^200 ) between the velocity dispersion in the 
dark matter distribution and the total mass within T2oo, we de- 
rive ax = V50 Ho r2oo = v50 Ho c r s for a NFW potential 
and compare our estimates of ax from the best-fit values in 
Table [I] to the optically-determined velocity dispersion, a v . 

Out of 16 clusters examined, we observe two systems 
(Al 19 and A754) lying with a relative difference in mass larger 
than 3 a. When the optically and X-ray determined veloc- 
ity dispersions are compared, three objects (A119, A754 and 
A2256) show significant deviations (see Fig. ||). These three 
clusters are NCF systems, are known to have irregular and 
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Fig. 3. Differences (in a) between the interpolated optical val- 
ues and the measured X-ray mass and velocity dispersion ver- 
sus the measured X-ray mass (see text for details). Filled circles 
represent CF galaxy clusters, whereas open squares are NCF 
objects. 

asymmetric X-ray brightness and, at least for A754 and A2256, 
are indeed involved in massive merging (e.g., A119: Ferretti 
et al. 1999; A754: Henriksen & Markevitch 1996; A2256: 
Molendi, De Grandi & Fusco-Femiano 2000) that may affect 
both the optical determinations of the velocity dispersion and 
the validity of the hydrostatic assumption made in the process 
of the estimation of the X-ray mass. 

In general, we measure a median deviation of about 1 .2 and 
0.8 a in mass and velocity dispersion measurements, respec- 
tively. Moreover, there is evidence that larger deviations are 
present in the subsample of NCF, not-relaxed systems (2.7 and 
2. Oct deviation in mass and velocity dispersion, respectively, 
for NCF; 1.0 and 0.7ct for CF). 

5. Relations among the observed quantities 

Our refined sample of 20 nearby (0.010 < z < 0.103; median 
redshift of 0.050) clusters of galaxies spans a factor of nore 
than three in mass-weighted temperature (3.1 keV < Tx < 9.9 
keV; median value: 6.9 keV) and two orders of magnitude in 
luminosity (1.7 x 10 44 erg s _1 < L x < 6.1 x 10 45 erg s _1 ; 
median value: 1.5 x 10 45 erg s _1 ). In the following analysis, 
we consider only the clusters with a total mass at a given radius 
larger than zero at the 95 per cent level of confidence, i.e. we 
select just the objects with a M /M < (1/1.96) = 0.51. 

In the present work, we investigate the correlations between 
the observed physical quantities in order to assess the robust- 
ness of the self-similar scaling relations for clusters of galaxies 
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Fig. 2. Comparison between deprojected (diamonds) observed values and best-fit /37-model (dashed line) for gas density, tem- 
perature and total gravitating mass. The dotted line indicates the best-fit temperature profile for a given mass model as described 
in Sect. 4. 



(e.g. Kaiser 1986). These relations are the product of simple 
assumptions on the formation and evolution of galaxy clusters. 
As non-linear structures, they are assumed to form by homoge- 
neous spherical collapse of gravitational instabilities of dark 
matter on which gas infalls are heated up by shocking pro- 
cesses. If no dissipation is considered, adiabatic X-ray emit- 
ting plasma can be considered to share the same potential well 
with dark matter with a spatial distribution that can be differ- 
ent from the dark matter's one but has to be the same at any 
earlier epoch (e.g. Bryan & Norman 1998, Arnaud & Evrard 
1999). However, deviations are expected from self-similarity 
under the effects of, for example, the dynamical history of clus- 
ters as three-dimensional aggregation of clumps (Jing & Suto 
2000, Thomas et al. 2001) and any additional physics acting 
on the intracluster gas over the simplistic infall in the poten- 
tial well (e.g. Evrard & Henry 1991; David, Forman & Jones 
1991; Bryan & Norman 1998; Bialek, Evrard & Mohr 2001; 
Borgani et al. 2002 and references therein). The latter case is 
particularly relevant to cool systems where the extra energetic 
amount required from their observed properties is comparable 
to their thermal energy (e.g. Ponman et al. 1996; Ponman, Can- 
non & Navarro 1999; Tozzi & Norman 2001). Our sample en- 
closes only clusters with temperature larger than about 3 keV 
and, thus, is expected not to be affected in a significant way 
by any increase of the gas entropy occurring during the clus- 
ter formation history. Therefore, we are able to investigate the 
galaxy cluster scaling laws excluding systematics related to the 
energetic budget. 



Table 3. Spearman's p rank correlation results on a set of phys- 
ical quantities. A small value in probability indicates signifi- 
cant correlation, na indicates the number of standard devia- 
tions by which the dependence in exam deviates from the null- 
hypothesis of uncorrelated data sets. The gas-mass-weighted 
temperature, T mw , is here used. 



relation 



Prob \na\ 
A =2500 



Prob I ncr I 
A =1000 



Ma-Ta 
Ra -T a 
La — Ta 
-Mgas,A — Ta 
/gas, a — Ta 
La - Ma 



0.69 
0.71 
0.65 
0.77 
0.58 
0.75 



0.001 
<0.001 

0.002 
<0.001 

0.007 
<0.001 



3.00 
3.09 
2.83 
3.34 
2.52 
3.26 



0.54 
0.54 
0.74 
0.79 
0.31 
0.64 



0.017 
0.017 
<0.001 
<0.001 
0.204 
0.003 



2.29 
2.29 
3.16 
3.35 
1.30 
2.70 



In our analysis, we adopt three different definitions for the 
plasma temperature at a given overdensity A: 

1 . the gas temperature in the shell at r A , 

T(r A ) = Ti(r A ), (8) 

2. the emission-weighted gas temperature within r A , 



-v0<r,<rA 



T ew (< r A ) = 



Li Ti 



■^0<ri<rA 



L(< r A ) 



(9) 
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3. the mass-weighted gas temperature within ta, 



T mw (< ?*a) 



E 



0<ri<rA 



M ■ T 

1 v± gas , l J i 



v-0<n<r A M rp 



E 



0<ri<r-A 



ill, 



gas : i 



AL 



gas (< r A ) 
(10) 

where i indicates the running cursor on shells and j on rings. 
We use the values of the temperature in the volume shells ob- 
tained from the best-fit procedure discussed in Sect. 4. In Fig. ^j, 
we show the temperature profiles defined above for a typical CF 
(A496) and NCF (A754) galaxy cluster. 

To check for the presence of linear dependence between the 
logarithmic values of the physical quantities being studied, we 
have performed both a non-parametric and a parametric analy- 
sis of our data. The two approaches can be considered as com- 
plementary: the non-parametric analysis has the advantage of 
not relaying on a specific model but does not treat errors; the 
parametric analysis, which assumes a power-law model, pro- 
vides a treatment for errors. 

For the non-parametric analysis we have used the Spear- 
man's p rank correlation of two sample populations (Press et 
al. 1992, p.634) . We quote in Table | the set of pair quan- 
tities with the respective Spearman's p and probability. Rela- 
tions with a ~ 3cr deviation from the null hypothesis of un- 
correlated datasets are present amongst the gravitating mass, 
the gas temperature, the gas luminosity and the gas mass. A 
weaker dependence appears between the gas mass fraction and 
the temperature. 

In our parametric analysis, we use the bisector modifica- 
tion (i.e. the best-fit results bisect those obtained from min- 
imization in vertical and horizontal directions) of the linear 
regression algorithm in Akritas & Bershady (1996 and refer- 
ences therein, hereafter BCES) that takes into account both any 
intrinsic scatter and errors on the two variables considered as 
symmetric. The uncertainties on the best-fit results are obtained 
from 10,000 bootstrap resamplings. 



5.7. M tl 



T relation 



In this and the following subsections we shall compare the 
the normalization and slope of the scaling relations obtained 
from our data to theoretical predictions based on the simplis- 
tic assumption of an isothermal sphere for both the gas rep- 
resented by its temperature and the collisionless dark matter 
particles (e.g. Kaiser 1986, Bryan & Norman 1998). As usual 
we indicate with 0t — (pm p a 2 ) / (fcT gas ), the ratio between 
the energy in the plasma and in the dark matter with veloc- 
ity dispersion, a. An isothermal distribution function is char- 
acterized by a proportional relation between the matter den- 
sity and the velocity dispersion, p{r) oc a 2 jr 2 (Binney & 
Tremaine 1987). This implies a total mass within a radius r of 
(2/3t &T)/(Gjum p ) r that can be compared to eqn. || to infer 
the relation between Af to t and T gas at given overdensity A: 



M tot (< r A ) 
lO^Me 



0.38 ,4 /2 



50 

IT 



1/2 



IkeV 



3/2 



logA/ 14 - -0.42 + 1.51ogT (+log ($ /a /n)) , (11) 



2.2 




Fig. 6. Plot of 1, 2, 3 a contour from x 2 statistic for the two 
interesting parameters of the linear fit applied to M — T mw 
(solid line), M - T cw (dashed line), M - T(R) (dotted line) 
relations at A =2500 and 1000 (the latter ones have higher 
normalization A). 



where H z is the Hubble constant at redshift z given in eqn. 0, 
fn = (^jf-^j (AS^) 1 / 2 and in the bottom row we have rewrit- 
ten the relation in log-log form. The assumption on the dark 
matter profile only affects the value of the normalization. We 
consider the isothermal case here adopted as a reference. 

In the considerations above, we have adopted the assump- 
tion that we are observing clusters just after their virialization 
(cfr. Voit & Donahue 1998 for the implication on the M — T re- 
lation of clusters that gradually form and stop evolving in a low 
density Universe). The M — T relation makes reasonable as- 
sumptions that have been tested both in numerical simulations 
and in observations. Moreover, this is a direct result coming 
from the combination of the conservation of energy through- 
out nearly-spherical collapse of clusters with the virial theorem 
(Afshordi & Cen 2002). 

In Table [| and Fig. |[ we show the results of the fitting 
analysis. A segregation is noticeable between CF and NCF ob- 
jects. When we fit the twelve CF clusters, we measure M14 = 
0.12(±0.06) x tAw 8(± ° 27) - When only the 8 NCF systems 
are considered, M u = 0.92(±2.01) x T°w 3(±100) . The slope 
does not show any significant change at the variation of the 
overdensity at which the quantities examined are considered 
(see Fig. 0). 

We compare now these results on the normalization and 
slope of the M— T relation with the values obtained in previous 
work. 
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Fig. 4. The temperature defined in eqn. 8, 9 and 10 are here plotted (dotted, dashed and dot-dashed, respectively) with the 
de-projected data points. 



Table 4. Results of the best-fit analysis. When the value of the slope is investigated, we apply the linear BCES bisector es- 
timator to the logarithmic of the power law Y = aX b , logF = A + BlogX (i.e. a = lO A ,b = B; errors in paren- 
theses). The temperature, T, is in unit of keV; the luminosity, L, in 10 44 /i^ 2 erg s 
the gas mass, M gas , in 10 13 /i7 5 M Q ; the radius at given overdensity, R, in 100 h 



50 



_1 ; the total mass, M, in 10 u h 5 oM Q ; 
kpc; the gas fraction, / gas , in h 5 ^ 2 . 



When the slope B is fixed, we estimate the median of the distribution of logy — BlogX. The scatter on Y is measured as 

2 1 1/2 

Y^j=\ n 0°S — A — BlogXj) /N . Note that the scatter along the X-axis can be estimated as cr log x — <?\ogY/B. 



relation 




A =2500 






^ =1000 










B 


o"io g y 


A 


B 


ClogY 


M14 


— T 

J- mw 


-0.70 (0.18) 


1.54 (0.22) 


0.14 


-0.52 (0.28) 


1.76 (0.34) 


0.25 






-0.60 (0.15) 


1.50 (fix) 


0.16 


-0.34 (0.20) 


1.50 (fix) 


0.24 


M14 


— T 

J ew 


-0.67 (0.20) 


1.51 (0.24) 


0.15 


-0.73 (0.27) 


1.94 (0.32) 


0.22 






-0.60 (0.16) 


1.50 (fix) 


0.16 


-0.39 (0.18) 


1.50 (fix) 


0.20 


M14 


-TOR) 


-0.61 (0.14) 
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Fig. 5. M — T mw relation. The solid and dashed lines represent our best-fit results for the given overdensity, using the slope 
as free parameter and fixing it to 1.5, respectively. Filled circles represent CF galaxy clusters, whereas open squares are NCF 
objects. (Left) The best-fit slope is 1.54(±0.22). The dotted line represents the best-fit result from Allen et al. (2001; A = 2500). 
(Right) Best-fit of the M — T cw relation at A = 500 (slope: 2.17 ± 0.37). The dotted line shows the best-fit from Finoguenov et 
al. (2001). 
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Fig. 7. (Left) Behaviour of the slope in the M — T mw , M gas — T mw , L — T mw and L — M relations as function of the considered 
overdensity. The diamonds + dashed error bars show the results for CF systems only. The dotted lines indicate the values predicted 
from the scaling laws. (Right) Values of the scatter of the same relations. Full dots represent the entire sample, whereas the 
diamonds indicate the results for CF objects. 
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In numerical simulations, fixing the slope to 3/2, the nor- 
malization Pt in eqn. [IT] ranges between 1.15 (model CL2 
in Navarro, Frenk & White 1995) to 1.24 (Evrard, Metzler 
& Navarro 1996) and - 1.3 (Bryan & Norman 1998; cfr. 
their Table 2). We measure a normalization (3t of 1.14(±0.34) 
that is consistent with the results quoted for simulated clus- 
ters. For example, we obtain a best-fit correlation of il/14 = 
0.46(±0.21) x (A = 1000), that has a normalization 
lower by 13 per cent (and only 0.3er apart) than the value mea- 
sured at the same overdensity in the gas-dynamic simulations 
ofEvrardetal. (1996). 

From an observational point of view, Horner, Mushotzky 
& Scharf (1999) claim that the M — T relation is steeper than 
the traditional scaling, following a oc 7 11 8-2 law, when the 
mass is estimated according to the j3— model. Slopes steeper 
than virial prediction are also observed in high-redshift clus- 
ters (Schindler 1999) and highly-luminous clusters (Ettori & 
Fabian 1999) samples, where isothermality is assumed. More- 
over, Neumann & Arnaud (1999) from a j3— model estimate of 
the gravitating mass obtain a M— T relation consistent with the 
classical scaling relation. Nevalainen, Markevitch & Forman 
(2000) from a sample of 6 clusters and 3 groups/galaxies with 
temperature profiles observed from ASCA and ROSAT , respec- 
tively, inferred at A = 1000 a slope of 1.79(±0.09) and a nor- 
malization significantly lower than the one observed in simula- 
tions, suggesting evidence of breaking of the self-similarity in 
the less massive systems due to heating processes. Finoguenov, 
Reiprich & Bohringer (2001) studied two samples of clusters, 
one comprising a complete sample of 63 bright objects from 
the ROSAT All Sky Survey with an assigned emission-weighted 
temperature collected from the literature and the other includ- 
ing 39 systems (22 of these with T < 3.5 keV) with known 
temperature profiles that are used to infer the total mass in 
combination with a /37-model (see Section 4.1 above). Cor- 
relating M to t(< r5oo) with an emission-weighted tempera- 
ture, these authors find a slope of 1.58(±0.07) for the sam- 
ple with a resolved temperature profile and excluding objects 
with A/500 < 5 x 1O 13 M (the slope increases slightly to 
1.78 ± 0.09 for the whole sample, consistent with the result 
from the flux limited sample). The normalization is more than 
50 per cent lower than the value quoted in Evrard et al. (1996). 
Allen, Schmidt & Fabian (2001), using spatially resolved X-ray 
spectroscopy with the Chandra observatory of five highly mas- 
sive cooling-flow (and so, relaxed) galaxy clusters at intermedi- 
ate redshifts, found consistency with the scaling law prediction 
and a normalization 40 per cent lower (but with a deviation sig- 
nificant only to 1.8(7 considering their quoted values) than what 
is observed in simulated clusters. 

It is worth noticing that a tight correlation between the pa- 
rameters A and B of the linear fit performed on the logarithmic 
values appears from the contour plot of the probability distri- 
bution obtained applying the \ 2 statistic (see Fig. ||). In this sit- 
uation, any slope larger than 1.5 (i.e. B > 1.5) requires a lower 
normalization A consistent with what is generally observed. 

We show in Fig. || a comparison with the recent results from 
Finoguenov et al. (2001) and Allen et al. (2001). Both are con- 
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Fig. 8. R — T mw relation. The solid line represents the best- 
fit slope of 0.47(±0.07); the dashed line is obtained fixing 
the slope to 1/2. Filled circles represent CF galaxy clusters, 
whereas open squares are NCF objects. 

sistent with our results. For example, Allen et al. measure a 
normalization that is 13 per cent lower than ours with a differ- 
ence of about 0.4cr. 

Another way to consider the relation between the gas tem- 
perature and a physical quantity related to the overdensity typ- 
ical of a cluster, is the R — T relation, where R is the radius 
of a sphere enclosing a given overdensity. It can be obtained 
directly from eqn. |ll] 

R _„i/2 /50^ ^1000^ 1/2 f T \ 1/2 



lOQ/ij^kpc 



5.07 py 2 



H z J \ A J 
0.5 log r(+ log (4 /2 /r 



(12) 



log -R100 = 0.71 -1- u.o iog 1 y-r rog yp T jn^ 

Fitting a power law as done above, we obtain i?ioo = 
3.27(±0.41) x T n °,w 7(±0 07) (A = 2500; cf. Table | and Fig. |). 
Fixing the slope to 1/2 in accordance with what predicted from 
scaling laws the normalization is of 530(±79) kpc (A = 1000), 
only 0.4cr below the value of 566 kpc measured in the hydro- 
dynamics simulations of Evrard et al. (1996). 



5.2. M gas — T relation 
The gas mass fraction 



/gasOA) 



A/gas (< r A ) 
Mtot(< r A )' 



(13) 



should be constant in a population of galaxy clusters that sat- 
isfies the self-similar behaviour. Once the self-similarity is 
broken, a dependence of / gas on the temperature is expected 




Fig. 10. (Left panel) M ga3 — T relation. The solid line represents the best-fit slope of 1.91(±0.29); the dashed line is obtained 
fixing the slope to 3/2. Filled circles represent CF galaxy clusters, whereas open squares are NCF objects. (Right panel) M gas — 
T C w relation. The solid line represents the best-fit M gaS) i3 = 0.52(±0.19) x Tew 89 ^ ' 20 "*; the dashed line is obtained fixing the 
slope to 3/2. The dotted line represents the best-fit results in Mohr et al. (1998). 



(Arnaud & Evrard 1999, Viklinhin et al. 1999). This propa- 
gates directly to the M gas — T relation in the following way: 
-Mgas ~ /gas Mtot oc T 3 ^ 2+a , where a dependence of the form 
/gas oc T a is introduced. In Fig. |9]and Table |], we show that 
there is a slightly positive correlation between / gas and T, with 
a slope larger than by 1.9a both at A = 2500 and 1000. As 
a consequence of this, the Af gas — T relation (Fig. |o|) tends to 
show a slope larger than 1.5. 

Fitting with a power law the outer X-ray emission, 
Vikhlinin, Forman & Jones (1999) found that Af gas oc 
T i.7i±o.i3 (at the baryon overdensity of 1000 » 200 in the dark 
matter overdensity). From the observed L — T relation and ap- 
plying X-ray scaling laws, Neumann & Arnaud (2001) found 
that Af ga s oc T 194 is required. Mohr, Mathiesen & Evrard 
(1999) estimated in a sample of nearby galaxy clusters that 
M gas oc T 198±018 within a density contrast of 500. At this 
overdensity, our sample includes 19 galaxy clusters, of which 
only two are directly observed (see Fig.|l]). If we extrapolate the 
physical quantities as described in Section 4, we obtain the cor- 
relation plotted in Fig. [l0| (right panel) that shows an agreement 
in the overall trend with a slope of 1.89(±0.20). 

It it worth noting that, while at an overdensity of 2500 there 
is a considerable segregation between CF and NCF systems, 
when going out to A =500 the segregation disappears. This 
kind of behaviour is in line with expectations indeed, while 
at A = 2500 the CF region contributes substantially to the 
overall gas mass, at A = 500 the contribution is quite small. 
The fact that segregation is actually not observed in our data at 



an overdensity of 500 indicates that the extrapolations made to 
derive quantities at A = 500 are reasonable ones. 

5.3. L — T relation 

The self-similar dependence between luminosity and tempera- 
ture is written as 

L w eVolK A(T) n 2 as R 3 

« / g 2 as T 1 ' 2 M 2 ot R- 3 « H 2 A f 2 as T 1 ! 2 M tot 

« H z A 1 / 2 / g 2 as T 2 , (14) 
where e is the X-ray emissivity and several, but reasonable, as- 
sumptions are made. Specifically: the cooling function, A(T), 
is here described by only bremsstrahlung emission that is 
strictly valid only for T > 2 keV; the gas fraction, / gas , and 
the radial dependence in the volume, Vol w r 3 ^, are not depen- 
dent on the temperature (see, e.g., Arnaud & Evrard 1999). 

On the other hand, observations of the X-ray properties in 
cluster samples did show a departure from the predicted slope 
of 2 with a measured value of about 3 (Mushotzky 1984, Edge 
& Stewart 1991, David et al. 1993) and a scatter along the 
mean relation that can be reduced considering properly the ef- 
fect of the cool cores (Fabian et al. 1994). Once the impact 
of the strength of cooling flows on the measured luminosity 
and emission-weighted temperature is taken into account, the 
L — T correlation is shown to be tighter and more, but still not 
completely, consistent with the self-similar prediction (Allen & 
Fabian 1998, Markevitch 1998, Arnaud & Evrard 1999, Ettori, 
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Fig. 9. /g as — T mw relation. The solid line represents the best- 
fit slope of 0.61(±0.31); the dashed line is obtained fixing the 
slope to 0. Filled circles represent CF galaxy clusters, whereas 
open squares are NCF objects. 

Allen & Fabian 2001). However, a steeper L — T dependence 
is expected when this relation is investigated over one order of 
magnitude in temperature, if the intracluster gas was pre-heated 
before the accretion in the potential well rising the entropy level 
in cooler systems (e.g. Ponman et al. 1996; Cavaliere, Menci & 
Tozzi 1997; Ponman, Cannon & Navarro 1999; Tozzi & Nor- 
man 2001). Gas-dynamic simulations (e.g. Bialek, Evrard & 
Mohr 2001, Borgani et al. 2002) are capable of recovering the 
observed slope once the effect of non-gravitational heat input 
from, e.g., AGNs and supernovae in the order of about 100 keV 
cm 2 , released either in an impulsive way or during the cluster 
formation history according to the star formation rate, is taken 
into account. 

In the temperature range investigated here, we observe a 
slope that is slightly higher, but still consistent within 2a, with 
2 using any of the three definitions of T (see Table |j) and at 
any overdensity (cf. Fig. ^|). 

5.4. L — M relation 

From the observed correlation between total mass and gas tem- 
perature and between X-ray luminosity and temperature, it is 
straightforward to derive the dependence between gas luminos- 
ity and total mass. Combining eqn. [I] and 

we obtain that 

L « H 2 A / g 2 as T 1 ' 2 M tot 

« A 7 / 6 f 2 as (15) 
Using a sample of 106 clusters observed with ROSAT 'PS PC 
in the energy range 0.5-2.0 keV and with total masses es- 
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Fig. 11. L — T mw relation. The solid line represents the best- 
fit slope of 2.79(±0.55); the dashed line is obtained fixing the 
slope to 2. Filled circles represent CF galaxy clusters, whereas 
open squares are NCF objects. The dotted line shows the best- 
fit from Allen et al. (2001), that follows the distribution of the 
CF clusters in our sample. 

A=2500 
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Fig. 12. L — M relation. The solid line represents the best- 
fit slope of 1.54(±0.26); the dashed line is obtained fixing 
the slope to 4/3. Filled circles represent CF galaxy clusters, 
whereas open squares are NCF objects. 
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timated through the /3-model, Reiprich & Bohringer (2002) 
measure a slope of 1.80(±0.08) in the Lboi — M200 rela- 
tion. This value is marginally consistent with our results en- 
closed between 1.84(±0.23), at A =2500, and 1.28(±0.15), 
at A =1000 (cf. Table |] and Fig. ^). These values lie between 
4/3, as predicted in the self-similar scenario, and 1 1/6, that is 
expected for a 'minimum-entropy' scenario in which cluster 
cores maintain signature of an earlier galaxy formation activity 
(cf. Evrard& Henry 1991). 

It is worth noticing that to use the X-ray galaxy clusters 
as cosmological probes and constraining the cluster mass func- 
tion, the L — T and M — T relations are generally applied to 
recover from the observed X-ray luminosity a corresponding 
mass. For a given luminosity, L, this procedure carries an in- 
trinsic scatter on the estimated mass, <Ti og Af> that is the convo- 
lution of the scatter on the temperature originated in the L — T 
relation, <Ti og T, with that present in the M — T relation. From 
our estimates, this intrinsic scatter is about 0.20 to be compared 
with the observed scatter in the L — M relation of 0.14, that is 
about 30 per cent less (we have considered results for emission- 
weighted temperatures at the overdensity of 2500; these values 
are 0.24 and 0.20 at A =1000, respectively). As a consequence 
of that, the direct application of the L — M relation as tight as 
observed provides a more robust way to infer the cluster masses 
from large X-ray surveys (see also discussion on the cluster 
mass function estimated using the L — M relation in Reiprich 
& Bohringer 2002). 

6. Summary and Conclusions 

We have used resolved gas temperature and density profiles 
obtained from direct deprojection of BeppoSAX spectral data 
to measure the gravitating mass profiles of 20 nearby (0.010 
< z < 0.103; median redshift of 0.050) clusters. This sample 
is, to date, the largest for which the physical quantities have all 
been obtained from the same dataset. 

At the overdensity (with respect to the critical density at a 
given redshift) of 2500, where 18 galaxy clusters in our sample 
have an observed temperature and luminosity, their gas mass 
weighted temperature spans over a factor of three in mass- 
weighted temperature (3.1 keV < Tx < 9.9 keV; median value: 
6.9 keV) and two orders of magnitude in luminosity (1.7 x 10 44 
erg s _1 < Lx < 6.1 x 10 45 erg s _1 ; median value: 1.5 x 10 45 
erg s -1 ). We recover the total gravitating mass adopting ei- 
ther a King or a NFW functional form of the potential well. 
We obtain that the X-ray mass estimates are generally in good 
agreement with optical measurements. After a selection that in- 
cludes in the analyses only systems with positive total mass at 
the 95 per cent level of confidence for a given overdensity A, 
we investigate in a consistent and robust way several correla- 
tions present between the quantities observed and derived from 
the X-ray analysis. Considering the high thermal energy asso- 
ciated with these objects, we are in the conditions to investigate 
the galaxy cluster scaling laws excluding systematics related to 
the energetic budget, but including statistical errors propagated 
from spatially-resolved gas density and temperature profiles. 



Considering two subsamples of clusters with (CF; 12 ob- 
jects) and without (NCF; 8 obj.) a mass deposition rate in the 
core (note that we do not consider any multi-phase gas in our 
spectral analysis, so that this classification only individuate sys- 
tems more -CF- or less -NCF- relaxed), we observe a segre- 
gation in the correlation planes according to this classification 
as originally suggested from Fabian et al. (1994) on the L — T 
relation. Moreover, the scatter measured in these relations is 
reduced when only CF systems are considered. These results 
point out an intrinsic scatter that is generally present in these 
correlations and is due to the dynamical status of the objects in 
exam. As expected, the scatter is greater for relations computed 
at larger overdensities, i.e. smaller radii, where the contribution 
from the more concentrated gas density is larger. 

For each scaling relation that we investigate in the present 
work, we summarize here our main results: 

- M — T relation: at any overdensity, the slope is consis- 
tent with a value of 1.5. In particular, we obtain the best-fit 
robust relations M 1A = 0.20(±0.08) x Tmw 4(±0 ' 22) and 
M14 = 0.25(±0.09) x T™. at A = 2500. A segrega- 
tion between CF and NCF clusters is observable. When 
the slope is fixed to 1.5, the scatter is reduced by about 30 
per cent for only CF objects (A = 2500). The normaliza- 
tion is lower than the results of hydrodynamic simulations, 
but the large scatter observed makes large the uncertain- 
ties that we measure. In the same direction, the results on 
the R — T relation show a slope in agreement with what 
is obtained in numerical simulations and a normalization 
slightly lower (R WQ = 3.27(±0.41) x T°w 7(±0 07) and 
i?ioo = 3. 19 (±0.32) x T,°-5). The evidence for a strong 
correlation between the fitted parameters in the M — T re- 
lation suggests that the physical interpretation of steeper 
slopes and lower normalizations can be misleading if the 
degeneracy between these quantities is not properly taken 
into account. 

- M gas - T relation: both at A = 2500 and 1000, we ob- 
serve a marginally significant trend between the gas mass 
fraction and the temperature (best-fit slope: 0.61 ± 0.31 at 
A =2500). As a consequence of that, there is evidence, 
significant at less than 2a, that the slope in the M gas — T 
relation is different (and more steep) than the one present 
in the M — T relation. At A = 2500, we measure a slope 
of 1.91 (±0.29). Due to the observed segregation between 
relaxed and not-relaxed clusters, the scatter is reduced by 
20 per cent when CF objects are considered. 

- L — T relation: our results are consistent with a slope 
of 2 as predicted from the scaling law relations (maxi- 
mal deviation of 2.\<j at A =1500). At the overdensity of 
2500, we measure L 44 = 0.06(±0.06) x T^ 9(±0 ' 55) and 
L 44 = 0.29(±0.22) x T n 2 lw ). The scatter is lowered by 20 
per cent when only CF galaxy clusters are considered. 

- L — M relation: the central value measured on the slope 
of this relation varies between the self-similar and the 
'minimum-entropy' prediction. The scatter in the relation 
is reduced by 25 per cent when only CF clusters are con- 
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sidered. We notice that the direct application of the L — M 
relation as tightly as observed provides a more robust way 
(in other words, less scatter is propagated) to infer the clus- 
ter masses from large X-ray surveys than the combined ap- 
plication of the L — T and M — T relations. 

In a forthcoming paper (Ettori et al., in preparation), the 
cosmological implications of these relations and of the gas 
mass fraction distribution observed at different overdensity will 
be investigated. 
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